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Based on the Konigl's inhomogeneous jet model, we estimate the jet parameters, 
such as bulk Lorentz factor F, viewing angle 9 and electron number density ric from 
radio VLBI and X-ray data for a sample of active galactic nuclei (AGNs) assuming 
that the X-rays are from the jet rather than the intracluster gas. The bulk kinetic 
power of jets is then calculated using the derived jet parameters. We find a strong 
correlation between the total luminosity of broad emission lines and the bulk kinetic 
power of the jets. This result supports the scenario that the accretion process are 
tightly linked with the radio jets, though how the disk and jet are coupled is not 
revealed by present correlation analysis. Moreover, we find a significant correlation 
between the bulk kinetic power and radio extended luminosity. This implies that the 
emission from the radio lobes are closely related with the energy flux transported 
through jets from the central part of AGNs. 
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1 INTRODUCTION 

The formation of highly relativistic jets in active galactic nuclei (AGNs) is one of the unsolved fundamental problems in 
astrophysics (e.g. Meier et al. 2001). It has been assumed that jets are produced close to the central black hole, involving 
power extraction from the black hole spin (Blandford & Znajek 1977; Macdonald & Thorne 1982; Thorne & Blandford 
1982) and/or from the accretion disk (Blandford & Payne 1982). Although the jet formation remains unclear, the estimate 
of the jet power is of fundamental physical interest, since it can be used to quantify the power emerging from the central 
engine of the radio source. Recently, the prescriptions for AGN feedback have been introduced into semi-analytic models of 
galaxy formation, and both Bower et al. (2006) and Croton et al. (2006) show that this feedback is able to solve the issue 
of the bright end of the luminosity function, whilst simultaneously solving other problems of galaxy formation models such 
as why the most massive galaxies are so red. Although the form of the AGN feedback adopted is very different in the two 
prescriptions, the relativistic ejecta from the AGN is a conceivably important ingredient of AGN feedback. Indeed, in clusters 
of galaxies containing powerful radio sources. X-ray observations have revealed bubbles and cavities in the hot intracluster 
medium, evacuated by the expanding radio source (e.g. McNamara et al. 2000; Fabian et al. 2003). Recent studies showed 
that the mechanical luminosity of radio sources are sufficient to suppress cluster cooling flows (Best et al. 2006; Nusser, Silk 
& Babul 2006). To understand the interaction between the radio sources and the surrounding medium, it is clearly important 
to estimate the bulk kinetic power of radio jets, since the expanding radio sources provides a direct way for the AGN output 
to be coupled to its environment. 

The relation between the jets and the accretion processes in active galactic nuclei has been extensively explored by many 
authors and in different ways. The strong correlations have been found between the low-frequency radio and narrow-line 
luminosities of 3C radio sources (Baum & Heckman 1989; Rawlings et al. 1989; Saunders et al. 1989), and also between the 
broad line and extended radio luminosity for radio-loud quasars (e.g. Cao & Jiang 2001). The link between the jets and 
the accretion processes can also been studied through exploring the relationship between luminosity in line emission and 
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kinetic power of jets in different scales (Rawlings & Saunders 1991; Celotti & Fabian 1993; Falcke, Malkan & Biermann 1995; 
Wang, Luo & Ho 2004). Rawlings & Saunders (1991) used the narrow-line luminosity as indicative of the accretion power and 
estimated the power transported by the jot from the energy content and lifetime of the radio lobes, finding a good correlation 
between the two. Using radio data on very long-baseline interferometry (VLBI) scales and the standard synchrotron self- 
Compton (SSC) theory, Celotti & Fabian (1993) estimated the jet kinetic power to put constraints on the matter content of 
jets. It offers some clues to understand the fundamental questions of the mechanisms, such as the coUimation and acceleration 
of jets. Celotti, Padovani & Ghisellini (1997, hereafter C97) explored the relation of luminosity in broad emission lines with 
the kinetic power of the jets for a sample of radio-loud AGNs. Their estimate of the bulk kinetic power is based on the 
adoption of the SSC model applied to the radio VLBI data and X-ray (or optical) fluxes. Lacking more accurate information, 
the minimum F for any given d [i.e. F = 0.5((5 -|- 1/(5)] is used in the derivation of bulk kinetic power for objects with 5 > 1, 
otherwise the F is derived from an average S. They found a suggestive hint of correlation between these two luminosities which 
is in favour of a link between the accretion process and the jets. However, by re-estimating the luminosity in broad emission 
lines on the sample of C97, Wang et al. (2004) argued that the jet bulk kinetic power is significantly correlated with the 
disk luminosity. Maraschi & Tavecchio (2003, hereafter MT03) found that the jet power is linearly proportional to the disk 
power for a sample of blazars, for which the jet powers were estimated using physical parameters determined from uniformly 
modeling their spectral energy distributions. However, by studying a sample of quasars from Wang et al. (2004), Punsly & 
Tingay (2005) argued that the bulk kinetic power and the bolometric luminosity are very weakly correlated in radio-loud 
quasars that possess blazar cores. 

In the framework of the relativistic beaming and the SSC model, the physical quantities in the jets can be estimated 
using the VLBI observations and the X-ray flux density. Marscher (1987) derived the beaming parameters on the assumption 
of homogeneous spherical emission plasma. Ghisellini et al. (1993) adopted Marschcr's approach and obtained the Doppler 
boosting factor 6 for 105 sources. Moreover, Readhead (1994) estimated the equipartition Doppler boosting factor 6eq, assuming 
that the sources are in equipartition between the energy of radiating particles and the magnetic field. Giiijosa and Daly (1996) 
derived the Scq for the same sample in Ghisellini et al. (1993). The variability Doppler factor 5var is derived on the assumption 
that the associated variability brightness temperature of total radio fiux density flares are caused by the relativistic jets 
(Lahteenmaki & Valtaoja 1999). The advantage of homogeneous sphere model is that the formalism is simple and the value 
of 5 derived is independent on the cosmology model. However, it is generally difficult to know the component angular size and 
the flux at the turnover frequency, so one has to assume that the VLBI observing frequency is the synchrotron self-absorption 
frequency. In addition, the dependence of core size on the observing frequency in some sources is inconsistent with the 
homogeneous spherical assumption. Blandford and Konigl (1979) and Konigl (1981) presented an inhomogeneous relativistic 
jet model, in which both the flat spectrum characteristics of some AGNs and the dependence of the core size on the observing 
frequency could be well explained. Based on their model, a new approach has been proposed to derive the jet parameters 
including bulk Lorentz factor F, viewing angle 9 and electron number density rie in the jets (Jiang, Cao & Hong 1998, hereafter 
J98). The proper motion measurements on the jets' components were adopted in their calculations. The correlation between 
the brightness temperature in the source rest frame and the derived Doppler factor suggested that the derived values of 
beaming parameters are qtiitc reliable (,J98). Moreover, the derived beaming parameters from the homogeneous sphere model 
is in general consistent with that from their inhomogeneous jet model. 

In this work, we foUow the method of J98 to derive the physical quantities of jets for a large sample of AGNs, then 
rc-analyze the relation between the luminosity in broad line emissions and the bulk kinetic power of the jets. In Section 2, we 
describe the sample of sources. The method of jet parameters derivation and the estimate of jet kinetic power are outlined 
in Section 3. Section 4 includes the results and discussion. In the last section, we draw our conclusions. The cosmology with 
Ho = 70 km s~^ Mpc~^, Q,m = 0.3, and Q,a = 0.7 have been adopted throughout the paper. 



2 THE SAMPLE 

In order to use inhomogeneous jet model to estimate the jet parameters, all sources should have VLBI measurements of 

proper motion of outflowing plasma. Combining with the relevant data of sources, such as the radio flux density, the size 
of the core and X-ray flux density, the jet parameters can be derived (J98). After searching the literature, our sample is 
constructed, which consists of 128 sources, including 94 quasars, 26 BL Lac objects and 8 radio galaxies. The observational 
data for the sample are presented in Table 1: (1) lAU name; (2) classiflcation of the source (Q= quasars; Qc= core-dominated 
quasars; Ql= lobe-dominated quasars; Qp= GHz peaked quasars; BL= BL Lac objects; G= radio galaxies); (3) redshift z; (4) 
observation frequency fs in GHz; (5) core radio flux density fc at frequency Vs; (6) VLBI core size 6d in mas; (7) reference 
for the VLBI data; (8) the proper motion /iapp; (9) reference for the proper motion; (10) 1 keV X-ray flux density /ikev in 
/iJy; (11) reference for the X-ray flux. 

When there are more than one moving components, we adopted the fastest one, which is regarded as a good approximation 
of jet bulk motion. In addition, we use the core flux density measured at the highest frequency, when VLBI core was measured 
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at more than one frequency. We assume that all the observed X-ray flux density is attributable to the SSC emission in the 
derivation, which will introduce some uncertainties. However, the derived jet parameters are not sensitive to the adopted X-ray 
flux density (J98). The rcdshift of 0716+714 is not available, and a value of 0.3 is assumed in the calculation. To calculate 
the total luminosity of broad emission lines, the available measurements of various broad emission lines for each source are 
collected from literatures. Moreover, we search the literatures and collect all available radio extended emission data from VLA 
observations for each source. The data are available for all 128 sources, of which the data from Australia Telescope Compact 
Array (ATCA) are used for the southern source 0208-512, and only upper limit is available for 8 sources due to the faintness 
or non-detection of extended emission. The extended flux density is k-corrected to 5 GHz in the rest frame of the source 
assuming ae = -1 (/ext <xi>"'). 



3 BULK KINETIC POWER 

We estimate the bulk kinetic power based on the inhomogeneous jet model (Konigl 1981). Using VLBI radio data, including 
proper motion of plasma, and X-ray fluxes, we can calculate the comoving electron number density rie, the magnetic fleld 
intensity, Lorentz factor T and the viewing angle 9. A brief description of the method is given below, and we refer to J98 and 
references therein for a complete description. In brief, the jet parameters were calculated by relating the model predicted size 
of optically think region in the jet, radio emission from the optically think region along the jet, SSC X-ray emission from the 
unresolved jet, and the apparent transverse velocity to the observables of radio core size, radio core flux density. X-ray flux 
density, and the proper motion. In this paper, we assume that the X-rays are from the jet rather than the intracluster gas 
and the radio blob speed is the jet flow speed. 

In Konigl's inhomogeneous jet model, the magnetic field B{r) and the number density of the relativistic electrons 
ne(r, 7c) in the jet are assumed to vary with the distance from the apex of the jet r as B{r) = Bi(r/ri)^"^ and ne{r,^e) = 
ni(r/ri)~"7e~*-^°''"^\ respectively, where ri = 1 pc and 7e is the Lorentz factor of the electron in the jet. Given that the bulk 
motion velocity of the jet is (3c (corresponding to a Lorentz factor F) with an opening half-angle 0, and the axis of the jet 
makes an angle 9 with the direction of the observer, the distance from the origin of the jet, r{Tvs = 1), at which the optical 
depth to the synchrotron self-absorption at the observing frequency Vs equals unity, is given as 

= (2c2(a)rim0csce)^/(^°+«)'=»(Bi5)(^«+^)/(^"+^"="'(i..(l + «))-^/'=» (1) 

ri 

where C2(a) is the constant in the synchrotron absorption coeflicient, 5 is the Doppler factor, and km = [2n -|- m{2a + 3) — 

2]/(2a-f 5). 

The projection of the optically thick region in the jet is then used as the observed VLBI core angular size 6d, 

e. = = (2) 

where 7?a is the angular diameter distance of the source. 

By integrating the omission from the optically thick region along the jet, the radio flux of the core can be obtained 

(3) 



_ r'Usme ci(q) p-i/2,,5/2 / ( rjr^, = 1) \ 



where i^s is the VLBI observing frequency, and ci (a) and C2 (a) are the constants in the synchrotron emission and absorption 
coefficients, respectively. 

Equation (13) in Konigl's work gives the X-ray fiux density estimation from an unresolved jet. As in J98, we adopt 
the expression in the frequency region Uc > Vcbirm), where ryi is the smallest radius from which optically thin synchrotron 
emission with spectral index a is observed (Konigl 1981). 

The proper motion observed with VLBI can be converted to the apparent transverse velocity /3app, which is related to 
the bulk velocity of the jet f3c and viewing angle 6, 

^''■'P ~ 1-/3 cose 

Given the three parameters a, rn, n, and the relation between the opening half angle (f) and the Lorentz factor F, 
the parameters of an inhomogeneous jet can be derived from VLBI and X-ray observations, using the above equations and 
equation (13) in Konigl (1981). In our calculation, we take a = 0.75, the opening half-angle 4> = l/F, and assume m = 1, 
n = 2 corresponding to a free jet (Hutter & Mufson 1986). 

With the estimated comoving total electron number density nt, Lorentz factor F and the cross section of the jet S, the 
bulk kinetic power is then derived as 

Lkin = 5nt(me(7) + m+(7+»c^F(F - l)/3c, (5) 
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where rUe is the electron rest mass, m+ is the rest mass of positive charge, (7) is the average Lorentz factor of electrons, and 
(7+) is the average Lorentz factor of positive charges. For a conical jet with an opening half-angle <f), S = 2nr'^{l — COSffl). 
The total electron number density nt is given by 

7max 

nt = 7le(r,7e)d7e, (6) 

Tmin 

The bulk kinetic power of the inhomogeneous jet becomes 
ikin = a-Vr-2-"rrni7-?„"(l - cos0)(me(7> +m+(7+))r(r - l)f3c^ (7) 
With our adoption of a = 0.75, (f> = l/F, m = 1, and n = 2, the bulk kinetic power of the jet is then given by 

Lkin = ^T^rlmj-lil - COS l/r)(me{7) + m+(7+))r(r - l)/3c^ (8) 

We note that the bulk kinetic power Lkin is independent of r, since n = 2 is adopted in the calculation and the particle 
conservation is then satisfied along r. Lkin is largely dependent of the matter content of jets and the low energy cut-off 7min 
of electrons. In present, the jet composition is still unclear, i.e. whether electron-positron or electron-proton (see Worrall & 
Birkinshaw 2006, for a recent review and reference therein). However, for an electron-proton plasma, 7min ~ 100 has been 
suggested, while 7min could be as low as unity for an electron-positron jet (e.g. Celotti & Fabian 1993). The detection of 
circular polarization strongly suggests that the jets arc electron-positron plasmas with 7min ^ 10 at least in some sources 
(e.g. Wardle et al. 1998). The similar conclusion is also arrived from powerful large scale X-ray jets, if they are interpreted as 
inverse-Compton scattering of cosmological microwave background photons in fast jets. From equation (8), however, we find 
that the bulk kinetic power Lkin for a electron-proton jet with 7min ~ 100 ('m+ = 1836 mo and (7+) = 1) is in agreement with 
that of electron-positron one with 7inin < 10 (m+ = me and (7+) = (7)) within a factor of three. In present work, we calculate 
the bulk kinetic power Lkin assuming electron-positron jets with 7niin = 1. A change to 7min = 10 will uniformly reduce Lkin 
by about a factor of three. However, in this work, we mainly focus on the correlation between the bulk kinetic power Lkin 
and the luminosity in broad emission lines Lblr, therefore, the value of 7min will not affect the correlation analysis, if the 
assumption that all sources have the same value 7min holds. 

The observational data necessary for calculations arc presented in Table 1. Following C97, we use the line ratios reported 
by Francis et al. (1991) and add the contribution from line H^ to derive the total broad line luminosity Lblr. The derived 
values of jet parameters, Lkin, Lblr and 5 GHz radio extended luminosity for our sources axe listed in Table 2: (1) lAU name; 
(2) the viewing angle of jet 6; (3) the Lorentz factor F; (4) the Doppler factor 5; (5) the normalization factor of electron 
energy distribution m; (6) the bulk kinetic power of jet Lkin; (7) the total luminosity in broad emission lines Lblr; (8) the 
references for flux of broad emission lines used to estimate Lblr; (9) the radio extended 5 GHz luminosity Lext.sGHz; (10) the 
references for the radio extended fiux. 



4 RESULTS AND DISCUSSION 

4.1 Bulk kinetic power and BLR luminosity 

Out of 128 sources, the measurements of various broad emission lines are only available for 98 sources from literature or 

Sloan Digital Sky Survey (SDSS) spectra, including 81 quasars, 15 BL Lac objects, and 2 radio galaxies. The reason of no 
measurements for remaining 30 sources could be either the non-detection of broad emission lines in 11 BL Lac objects and 
6 radio galaxies, or no published line flux measurements for 13 quasars. We have calculated the total luminosity of broad 
emission lines for these 98 sources. The relationship between Lblr and Lkin is shown in Fig. 1. We find a strong correlation 
between these two luminosities with a Spearman correlation coefficient of r = 0.565 at ^ 99.99% confidence. It should be 
noted that this correlation may be caused by the common dependence on redshift. We present the bulk kinetic power and 
BLR luminosity as functions of redshift z for the sample in Fig. 2. We therefore use the partial Spearman rank correlation 
method (Macklin 1982) to check this correlation. Still, a significant correlation with a correlation coefficient of 0.323 is present 
at about 99.9% significance level between Lblr and Lkin, independent of the redshift. We also perform a statistic analysis 
on the sources in the restricted redshift range 0.5 < z < 1.0. For this subsample of sources, we check the correlation between 
luminosity and redshift, and no correlation between either the bulk kinetic power or the BLR luminosity and redshift, is found 
(see Fig. 2), while a significant correlation is still present at 99.1 per cent confidence between the bulk kinetic power and total 
broad-Iinc luminosity (sec Fig. 1). Therefore, we conclude that this correlation might be intrinsic, at least for our present 
sample. Assuming that the BLR luminosity is due to the reprocessing of the ionizing radiation from the accretion disk, it 
therefore strongly supports the scenario of a tight connection between the relativistic jet and the accretion process. 
For all 98 sources, the ordinary least-squares (OLS) bisector method gives the following fit in Fig. 1: 
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log ikin = (0.86 ± 0.07) log iBLR + (8.78 ± 3.05) (9) 

Apart from finding a significant correlation, Rawlings & Saunders (f 991) found that the relationship between the bulk kinetic 
power and narrow line luminosity is close to proportionality, Q oc L^^^"'^, which extends over four orders of magnitude. 
Our relation of Lkin oc iBLR*" consistent with the relationship between the jot power and narrow-line luminosity in 
Rawlings & Saunders (1991) and Celotti & Fabian (1993), whereas it is somehow deviated, but not much, from the linear 
relation between the jet power and disk power found by MT03. However, it is much steeper than that of Wang et al. (2004), 
Lkin OC Lblr- Although what cause these differences is unclear, we note that the methods used to estimate the jet power 
in MT03 and Wang et al. (2004) are different from ours. The jet powers of MT03 were estimated using physical parameters 
determined from uniformly modeling their spectral energy distributions, while Wang et al. (2004) directly used the jet bulk 
kinetic power from C97, which was estimated using the homogeneous sphere SSC model. Moreover, MT03 obtained their disk 
luminosities either directly from the optical-UV luminosity of the big blue bump or from the original prescription of C97. The 
method Wang et al. used to estimate Lblr is basically same as ours. In addition, Wang et al. sample consists of 35 blazars, 
and only 16 sources (11 quasars and 5 BL Lac objects) were considered in MT03. Their samples are much smaller than our 
sample. Whether these factors influence the dependence of Lkin on Lblr needs further investigations. Despite this, our results 
strongly support the scenario that the accretion process are tightly linked with the kinetic power in the jet, though how the 
disk and jet arc coupled is not revealed by present correlation analysis. 

In general, BL Lac objects, which is thought to be FR I radio galaxies pointing at us, are characterized by very weak 
or absent emission lines, invisible blue bumps, and relatively powerful jets. Prom Fig. 1, it is clear that BL Lac objects 
have fainter broad-line luminosity compared to quasars, though only 15 BL Lac objects arc in our sample. We find that the 
Lkin — Lblr relation of BL Lac objects deviate from that of quasars, although it generally follows that of the whole sample. 
The linear fit using OLS bisector method for BL Lac objects shows 

log Lkin = (0.71 ± 0.12) log Lblr + (15.85 ± 4.99) (10) 
while for quasars, we have 

log Lkin = (1.12 ± 0.08) log Lblr - (2.94 ± 3.39) (11) 

Although the mechanism of jet formation is unclear, the different dependence of Lkin on Lblr in BL Lac objects and quasars 
can be due to the difference of the accretion power as measured in units of the Eddington one. Compared to quasars, BL 
Lac objects are characterized by radiatively inefficient accretion disks (Cao 2003), thus in these sources the jet power may be 
relatively dominant. To further check the Lkin — Lblr correlation in Fig. 1, we re-examine it for quasars only. When BL Lac 
objects and radio galaxies are excluded, we still find a strong correlation between Lblr and Lkin with correlation coefficient 
of r = 0.380 at 99.95% confidence. This further confirms the tight link between the accretion process and the kinetic power 
in the jet. 

If the accretion process and the jet formation are indeed closely related, then the tight relation between the mass 
channelled into jets and that accreted by black hole would be expected. We then investigate the relationship between the 
mass outflowing rate and accretion one, on assumption of Lboi IOLblr (Netzer 1990), and with the expression of the kinetic 

and the accretion powers as 

Lkin = FAfoutc^ Lbol = ^^MinC^ (12) 

where Mout is the mass outflowing rate, F is jet Lorentz factor, Min is the mass accretion rate, and rj is the efficiency of 
mass to energy conversion for accretion. Adopting the typical value of rj ~ 0.1, we find a significant correlation between Mout 
and Min for whole sample, with a Spearman correlation coefficient 0.514 at ^ 99.99% confidence. This implies that the mass 
outflowing rate in jet is closely linked with the accretion one in accretion disk. 

The present analysis is based on the derivation of jet parameters using inhomogeneous jet model. Some parameters and 
assumptions are adopted in the inhomogeneous jet model to derive the physical quantities of the jet (J98), which may induce 
some uncertainties in the estimation of Lkin. The most important is probably the intrinsic differences of the low energy 
cut-off 7min of clcctrons between the radio sources themselves. In this work, we adopt the same value of 7niin in deriving the 
kinetic power of the jet, which may not be true. However, we are not able to estimate 7min for each source at present stage. 
Nevertheless, we believe that the adoption of an uniform 7niin for all sources in the correlation analysis would not affect the 
main conclusion drawn here. Moreover, our results are based on the assumption of a = 0.75, m = 1, and n = 2, and we 
adopted these same values for all sources in our model calculation. However, we flnd that the alternative adoption of a, m, 
and n do not change our main conclusion, e.g. the strong correlation between Lkin and Lblr. In practice, the sources may 
have different values of parameters a, m, and n, and, in principle, these parameters could be constrained by the observable 
quantities (,198). Unfortunately, the information is only found for a few cases through multi- frequencies VLBI observations. 
Further high resolution multi-frequencies VLBI observations would be helpful to improve our model calculations. 

It should be noted that not all sources in our sample have available BLR luminosity. Therefore, the selection effects may 
be introduced in our correlation analysis, i.e. those sources without published broad line fiux measurements may likely be 
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biased towards those with weak lines, especially the 11 BL Lac objects without broad line flux measurements. To evaluate 
the selection effects, we tentatively calculate the upper limit of BLR luminosity assuming equivalent width of broad H/3 
lino EW < 5A for these 11 BL Lacs. Moreover, wo calculate the BLR luminosity for 13 quasar without BLR luminosity by 
randomly assigning the BLR flux in the BLR flux range of 81 quasars with BLR luminosity. Combining these 24 sources with 
those having BLR luminosity, we use the Astronomy Survival Analysis (ASURV) package (Isobe, Feigelson & Nelson 1986) 
to investigate the correlation by taking the upper limit into account. A significant correlation is still found with Spearman's 
rho correlation method. This correlation is confirmed by using the partial correlation method for censored data of Akritas 
& Siebert (1996) to exclude the common dependence of redshift. Furthermore, the correlation remains significant even we 
conservatively adopt broad H/3 line EW < lA for 11 BL Lacs. It thus seems that the non-BLR luminosity sources do not 
affect our correlation results. 

4.2 Bulk kinetic power and radio extended luminosity 

In Fig. 3, the relation between bulk kinetic power and 5 GHz radio extended luminosity is shown for all 128 sources. Since 

only the upper limit of 5 GHz extended luminosity is given for 8 sources, wc use ASURV package (Isobe ct al. 1986) to 
investigate the correlation and perform the linear regression analysis for our censored data. We flnd the significant correlation 
with a Spearman's rho correlation coefficient of r = 0.493 at » 99.99% confidence. Using the partial correlation method for 
censored data of Akritas & Siebert (1996) to exclude the common dependence of redshift, a significant correlation is still 
present between Lext.sGHz and Lkin. We use the Schmitt-binning method (Schmitt 1985) to perform y/x and x/y fits and 
then calculate a bisector of these two fits, as described in Shapley, Fabbiano & Eskridge (2001) (see also Isobe et al. 1990). 
We obtain: 

log Lkin = (0.82 ± 0.09) log Lext,6GHz + (12.38 ± 4.00) (13) 

which is shown as the solid line in Fig. 3. 

The extended radio flux is usually emerged from the optically thin radio lobes, and thus is free from the Doppler boosting 
effects, since the lobe material is generally thought to be of low enough bulk velocity. Therefore, the extended radio luminosity 
can be a good tracer of jet power (e.g. Cao & Jiang 2001). The significant correlation between Lkin and Lcxt.sGHz implies that 
the emission from radio lobes are tightly related with the energy ejected into the jet from the central parts of AGNs. This 
result is not surprising, as it could be naturally expected. Although the detailed mechanism of jet formation is still unclear, 
the energy can be transported through the jets to the radio lobes once the jets are generated. Most of the energy fiux of jets 
is not radiated away, instead are in mechanical form (i.e. bulk kinetic power), of which the particles and fields are necessary 
to produce the synchrotron luminosity that is detected in the radio lobes. If the radiative efficiency of radio lobes are similar 
between our radio sources, then the tight link between the bulk kinetic power of jets and the radio extended emission is 
expected, since the latter is optically thin and not effected from the Doppler enhancement. 

Motivated largely by the observed effects of radio-loud AGN on their environments at galaxy cluster scales (e.g. Fabian et 
al. 2003) , whether the heating effect of AGN activity, particularly radio-loud AGN activity, can balance the cooling of the gas 
has recently arose much interest (e.g. Best et al. 2005, 2006; Croton et al. 2006; Nusser, Silk & Babul 2006). Independently of 
the radio properties, Birzan et al. (2004) estimated the mechanical luminosity associated with the radio source, by studying 
the cavities and bubbles that are produced in clusters and groups of galaxies due to the interactions between the radio sources 
and the surrounding hot gas. The dependence of the mechanical luminosity on the 1.4 GHz radio luminosity of the associated 
radio sources, Lmeoh oc L°^^^° °^ fitted for their entire sample, is somewhat deviated from ours (equation 13). However, the 
dependence for the radio-filled cavities only (see Birzan et al. 2004 for details), Lmoch oc L'^i^^^^'^ is marginally consistent with 
ours within the errors. Despite this, we note that their work is mainly based on the galaxy clusters and the radio sources in 
these clusters, however, our present study focused on the powerful radio sources. Moreover, the jet power is estimated using 
different methods, and the different radio luminosity is used. In present, we are not able to draw a solid calibration between 
the radio emission and the kinetic power of jets, and it needs further investigations. 

Despite the strong correlation presented in Fig. 3, the significant scatter is clearly seen. This is not surprising since even 
for a source of fixed jet kinetic power the radio luminosity changes as the source ages (e.g. Kaiser et al. 1997). However, there 
are several factors that can introduce the scatter into the correlation. We note that the observed radio extended emission has 
been dissipated over a long period, which is not contemporaneous with the estimated bulk kinetic power. Moreover, when jets 
transported the energy flux from the central parts of AGNs to outer radio lobes, the jets can be decelerated by the interaction 
with the nuclear ISM and/or the entrainment of external gas (e.g. Tavecchio et al. 2006). As a result, the part of jet power 
will be lost to the ISM. In some extreme cases, the kinetic power of the jet on kiloparsec scales could be about three orders 
of magnitude weaker than the power of the jet on 10 - 100 pc scales duo to the jct-ISM interaction, i.e. virtually all of the 
jet power can be lost to the ISM within the inner kiloparsec (Gallimore et al. 2006). Consequently, the difference of jet-ISM 
interactions between the radio themselves may bring scatter. Furthermore, part of the scatter may be due to the different 
radiative efficiency in individual source. 
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4.3 Bulk kinetic power versus radiative luminosity 

It is well known that the monochromatic radio luminosity does not provide a good indicator of the mechanical energy output 
of a radio source. Radio sources are inefficient radiators. Bicknell (1995) estimated that the kinetic energy output of a radio 
jet is typically a factor of 100-1000 higher than the total radio luminosity of a radio source, which is recently confirmed by 
the observations of Birzan et al. (2004). 

We estimate the amount of radiative dissipation on parsec scales, i.e. the ratio of Z/kin/irad,in of the bulk kinetic power 
to the intrinsic radiative luminosity. The latter has been computed from the observed VLBI radio core fluxes. The Doppler 
correction on the monochromatic luminosity is assumed to be Lobs = S'^ Lint, where Lobs and Lint are the observed and intrinsic 
(comoving) luminosities. We can calculate p from the dependence of the core radio flux on the Doppler factor from equation 
(3), in which p is dependent of the value of a, m and n. In Fig. 4, we show the histogram of the derived ratios Lkin/Lrad.in- 
We find that the ratio covers about three orders of magnitude with the average value < log (Lkin/Lrad,in) >= 4.98 ± 0.79. 
This result is consistent with that of Celotti & Fabian (1993), although a wider spread of ratio distribution in their sources. 
The results indicate that for all sources the kinetic power is dominant with respect to the radiative output, and consequently 
that the radiative dissipation is not an efficient process. Moreover, the large variation in this ratio indicates that the radio 
luminosity is not necessarily a reliable probe of the available bulk kinetic power. 

It is commonly accepted that the synchrotron emission of blazars can extend to optical and even X-ray region, which 
can dominate over the thermal emission from accretion disk, and the radio emission as well. The spectral energy distribution 
(SED) of blazars is usually composed of two peaks, of which the first one is dedicated to the synchrotron emission for jets, 
and the second is due to the inverse Compton process (Fossati et al. 1998; Ghisellini et al. 1998). Thus the total radiative 
luminosity of the jets, if we can estimate from integrating over the synchrotron domain of SED, can represent the minimum 
power that must be associated with the jet in order to produce the observed luminosity. In this sense, the radio emission 
solely might not be a good indicator of the radiative output from the radio jets, i.e. the radiation losses of the kinematic jet 
flow. MT03 found that the radiative efficiency, i.e. the ratio of the total radiative luminosity of the jet to the jet power, can 
be in the range 1%-10%. Even so, the radiative dissipation is still not an efficient process, and the most of the energy flux is 
in the kinetic form. In present, it is not readily to estimate the total radiative luminosity of jets for our sample. Nevertheless, 
we believe that the inefficient radiators of radio sources would be still retained. 



5 CONCLUSIONS 

Based on the inhomogeneous jet model, we have calculated the jet parameters for a sample of AGNs. The bulk kinetic power 
of radio jets are then estimated using the derived jet parameters. We found a significant correlation between the bulk kinetic 
power of the relativistic jet and the total luminosity in broad emission lines, implying a tight link between the jet and accretion 
process. Moreover, the bulk kinetic power of jets are strongly correlated with the radio extended luminosity. This indicates a 
closely connection between the emission from radio lobes and the energy flux transported through jets from the central parts 
of AGNs. In addition, we found that the bulk kinetic power is dominant with respect to the radiative output, which means 
the radiative dissipation is not an efficient process. 
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Figure 1. The bulk kinetic power versus BLR luminosity. The circles represent quasars, and the triangles are BL Lac objects, while the 
rectangles show radio galaxies. The filled symbols arc the sources in the redshift range 0.5 < z < 1.0. The solid line is the fitted line for 
the whole sample using the OLS bisector method, and the dashed line is fitted for quasars only, while the dot-dashed line is fitted for 
BL La<; objects only. 
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Figure 2. The bulk kinetic power versus redshift (upper panel) and the BLR luminosity versus redshift planes (lower panel) for the 
sample. The circles represent quasars, and the triangles are BL hax: objects, while the rectangles show radio galaxies. The restricted 
redshift range, 0.5 < z < 1.0, is indicated with the dotted lines. 
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Figure 3. The bulk kinetic power versus radio 5 GHz extended luminosity. The symbols are the same as in Fig. 2. The arrows indicate 
the upper limit of extended luminosity. The solid line is the bisector linear fit using Schmitt-binning method for censored data (see text 
for details). 




Figure 4. The histogram shows the distribution of the ratio of the bulk kinetic power to intrinsic core radio radiative luminosities. 
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Table 1. VLBl and X-ray Data of tlit- Sample. 
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Table 1. Continued. 
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Table 1. Continued. 



Source 


Type 


z 
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11 
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Column (1): lAU name; Column (2): classification of the source (Q= quasars; Qc= core-dominated quasars; 
Ql= lobe-dominated quasars; Qp= GHz peaked quasars; BL= BL Lac objects; G= radio galaxies); Column 
(3): redshift z; Column (4): observation frequency i^s in GHz; Column (5): core radio flux density /c at 
frequency Ug; Column (6): VLBI core size 9^ in mas; Column (7): reference for the VLBI data; Column (8): 
the proper motion /lapp; Column (9): reference for the proper motion; Column (10): 1 keV X-ray flux density 
/ikcV ill /^Jy; Column (11): reference for the X-ray flux. 

References: (1) Kovalcv ct al. (2005); (2) Brunthalcr ct al. (2000); (3) Fey & Chariot (1997); (4) J98; (5) 
Jorstad et al. (2001); (6) Rokaki ct al. (2003); (7) Lister (2001); (8) Fey & Chariot (2000); (9) Hough ct al. 
(2002); (10) Fey, Clcgg & Fomalont (1996); (11) KcUcrmann ct al. (2004); (12) Vcrmculcn & Cohen (1994); 
(13) Siebert ct al. (1998); (14) Brinkmann, Yuan & Sicbcrt (1997); (15) Donato et al. (2001); (16) Ghisellini 
et al. (1993); (17) Comastri ct al. (1997); (18) Gambill ct al. (2003); (19) Marshall et al. (2005); (20) Reich 
et al. (2000); (21) Imanishi & Ueno (1999); (22) Galbiati et al. (2005); (23) Bloom et al. (1999); (24) Donato, 
Sambruna & Gliozzi (2005). 
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07'^4-i-i nn 


6.6 


12.8 


8.1 




46.38 






40.39 


30 


0J^04-I-4QQ 


6.9 


8.8 


8.3 


R R7P-I-04 


47.24 


45.39 


1^ 


42.04 


25 




1.9 


9.6 


17.3 


l.llE-i-04 


46.46 






43.59 


30 


nsi 4-1-42 


5.4 


6.2 


9.2 


111 p-uns 


45.43 


41.81 




41.34 


25 




4.2 


14.6 


13.7 




46.36 


43.40 


1 


40.82 


25 


0897-1-94^ 


4.0 


32.1 


10.7 


1 9'^P-I-O^ 


47.54 


44.93 


g 


42.18 


32 


089Q-I-04R 

u ii c »-' ^ "J 


4.8 


17.6 


11.3 


1 Q7P_l_n'R 


45.74 






40.97 


30 


08^R-I-71 
uooun^ / xu 


4.8 


29.3 


8.5 


1 PiOP-l-OR 


48.63 


46.43 


1^ 


43.78 


25 


08'^04-'S81 


8.7 


37.6 


2.2 


9 09P-knR 


48.76 


45.66 




43.78 


25 


0851+202 


6.3 


14.3 


8.3 


1.87E+04 


46.71 


43.60 


1 


40.48 


25 


0859-140 


3.9 


16.8 


14.7 


4.13E+04 


47.06 


45.74 


1 


43.45 


25 


0906+015 


6.3 


12.2 


8.7 


7.05E+04 


47.28 


45.11 


1 


42.53 


25 


0906+430 


1.8 


11.0 


19.5 


2.95E-^03 


45.90 


43.35 


1 


43.86 


25 


0917+449 


3.7 


13.1 


15.2 


8.70E+04 


47.37 


45.21'* 


9 


44.66 


33 


0917+624 


2.5 


16.5 


21.8 


1.91E+04 


46.72 


45.06 


1 


42.14 


25 


0923+392 


15.9 


43.4 


0.6 


1.77E-^07 


49.71 


45.78 


1 


43.45 


25 


0945+408 


2.4 


22.7 


23.8 


1.55E+04 


46.64 


45.59 


1 


43.15 


25 


0953+254 


8.4 


22.4 


3.8 


1.39E+05 


47.59 


44.97 


1 


41.85 


25 


0954+658 


6.8 


10.3 


8.3 


4.17E+03 


46.04 


42.63 


1 


41.32 


25 


1012+232 


4.9 


9.3 


11.3 


5.70E+03 


46.17 


45.16 


10 


42.73 


34 


1015+359 


7.6 


15.2 


6.1 


1.07E+05 


47.47 


45.98 


9 


<43.48 


35 


1039+811 


7.4 


4.9 


7.0 


2.53E-^04 


46.76 






42.67 


23 
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Table 2. Continued. 



Source 

1 




(degree) 




r 




S 




ni 
(cm"" 


') 


log 
(erg 


'^kin 


log -Lblr 
(erg s~i) 


Refs. 


log Loxt,,5GHz 
/ — 1 \ 

(erg s ^) 


Refs. 


1040+123 


18, 


.0 


11, 


.8 


1, 


.6 


1, 


.lOE+06 


48. 


.47 


45, 


.11 


1 


44.04 


25 


1049+215 


9 


.5 


10, 


.0 


5, 


.3 


1, 


.35E- 


^05 


47. 


,55 








44.00 


36 


1055+018 


1, 


.9 


5 


.5 


10, 


,7 


2, 


.80E- 


^04 


46. 


,82 


44, 


.53 


1 


43.13 


25 


1055+201 


10, 


.0 


15, 


.9 


3, 


,7 


2 


.19E- 


^05 


47. 


78 


46, 


.09 


11 


44.31 


30 


1101+384 


29, 


.3 


3 


.3 


1, 


,8 


1, 


.67E- 


^03 


45. 


,51 


41, 


.40 


12 


39.78 


28 


1127-145 


3, 


.5 


69, 


,7 


7, 


,2 


9, 


.81E- 


^05 


48. 


,45 


45, 


.77 


1 


43.21 


25 


1128+385 


6, 


.6 


2, 


.1 


3, 


.7 


7, 


.96E+04 


47. 


,01 


46, 


.26 


9 


42.26 


30 


1137+660 


45, 


.7 


5, 


.1 


0, 


.6 


6, 


.98E+05 


48. 


,21 


45, 


.83 


1 


43.94 


25 


1150+812 


15, 


.1 


11. 


.2 


2. 


.3 


6, 


.81E+05 


48. 


.26 








42.57 


30 


1156+295 


2, 


.5 


22, 


.2 


23. 


,1 


1, 


.05E+04 


46. 


.47 


44, 


.90 


5 


42.99 


5 


1219+285 


25, 


.1 


7, 


.0 


1, 


,4 


1, 


.12E- 


^04 


46. 


,44 


42, 


.25 


13 


38.94 


28 


1222+216 


3, 


.6 


27, 


,4 


13, 


,9 


1, 


.21E- 


^04 


46. 


,54 


44, 


.73 


14 


43.01 


21 


1226+023 


6, 


.2 


27, 


.3 


5, 


,6 


1, 


.29E- 


^05 


47. 


,56 


45, 


.59 


1 


43.14 


25 


1228+127 


85, 


.3 


1, 


.6 


0, 


,7 


1, 


.21E- 


^02 


44. 


,00 








40.84 


23 


1253-055 


4, 


.1 


10, 


.4 


13, 


,5 


3, 


.83E- 


1-04 


47. 


,01 


44, 


.64 


1 


43.46 


25 


1302-102 


14, 


.3 


6, 


.3 


3, 


.7 


1, 


.33E+04 


46. 


,51 


44, 


.91 


1 


41.83 


25 


1308+326 


2, 


.2 


45, 


.3 


22. 


.5 


4, 


.41E+04 


47. 


,10 


45, 


.12 


1 


42.56 


25 


1334-127 





.5 


9 


.8 


19. 


.4 


1, 


.06E+04 


46. 


,45 


44, 


.18 


1 


42.67 


25 


1345+125 


75, 


.3 


6, 


.4 


0. 


.2 


5, 


.13E+05 


48. 


,10 








<41.71 


35 


1406-076 


2 


.2 


61, 


,1 


18. 


1 


1, 


.OlE- 


^05 


47. 


,46 


45, 


.47 


15 


<43.81 


33 


1413+135 


2, 


.1 


10, 


,7 


18. 


,6 


5, 


.48E- 


^02 


45. 


,16 








<39.06 


37 


1458+718 


6, 


.4 


13, 


,1 


8, 


,3 


4, 


.58E- 


^04 


47. 


,09 


45, 


.47 


1 


44.74 


25 


1508-055 


3, 


.1 


43, 


,6 


13, 


,4 


8, 


.48E- 


1-04 


47. 


,39 


45, 


.52 


15 


44.72 


36 


1510-089 


5, 


.2 


27, 


,9 


7, 


,5 


4, 


.26E- 


f04 


47. 


,08 


44, 


.64 


1 


42.25 


25 


1532+016 


8, 


.0 


49, 


.7 


2, 


.0 


5, 


.15E+06 


49. 


,17 


44, 


.84 


1 


42.51 


25 


1546+027 


7, 


.7 


2, 


.4 


4. 


.2 


1, 


.85E+04 


46. 


.46 


44, 


.68 


1 


41.27 


25 


1548+056 


14 


.3 


3, 


.6 


4. 


.0 


5, 


.14E+04 


47. 


.02 








44.27 


36 


1606+106 


3 


.2 


42, 


.8 


12. 


.5 


1, 


.57E+05 


47. 


.65 








42.58 


30 


1611+343 


2 


.9 


82, 


,9 


9, 


,1 


1 


.07E- 


^06 


48. 


,49 


45 


.91 


1 


43.39 


25 


1618+177 


33, 


.1 


10, 


,3 


0, 


.6 


1, 


.17E- 


^06 


48. 


,49 


46, 


.14 


5 


43.46 


5 


1622-297 


4, 


.9 


21, 


,3 


9, 


,8 


6, 


.06E- 


^04 


47. 


,23 








<43.43 


35 


1633+382 


4, 


.2 


15, 


,7 


13, 


,7 


9, 


.90E- 


1-04 


47. 


,43 


45, 


.84 


1 


43.04 


25 


1637+826 


77, 


.7 


1, 


.1 


1, 


,0 


4, 


.94E- 


f02 


43. 


,76 








40.04 


23 


1641+399 


6, 


.1 


29, 


.5 


5, 


.4 


5, 


.91E+05 


48. 


,23 


45, 


.47 


1 


43.05 


25 


1642+690 


2, 


.5 


16, 


.7 


21. 


.7 


6, 


.14E+03 


46. 


.23 


43, 


.86 


1 


42.98 


25 


1652+398 


45, 


.0 


4, 


.2 


0. 


.8 


1, 


.24E+04 


46. 


.43 


41, 


.36 


12 


39.43 


28 


1655+077 


6, 


.8 


27. 


.7 


4. 


.7 


2, 


.28E+05 


47. 


.81 


43, 


.62 


5 


42.76 


5 


1656+053 


24, 


.2 


7, 


.9 


1. 


.3 


1, 


.12E+06 


48. 


.45 


46, 


.26 


5 


42.78 


5 


1656+477 


11 


.1 


4, 


.5 


5. 


.1 


5, 


.54E+04 


47. 


.09 


45, 


.76 


16 


<43.51 


33 


1721+343 


29, 


.7 


14, 


.5 


0. 


.5 


7, 


.05E- 


1-05 


48. 


.28 


44, 


.62 


1 


42.22 


25 


1730-130 


3, 


.6 


26, 


,2 


14. 


1 


1, 


.32E- 


^05 


47. 


.57 


44, 


.83 


17 


42.93 


32 


1749+096 


1, 


.1 


8, 


.9 


17. 


,3 


2, 


.62E- 


^03 


45. 


,83 








39.74 


28 


1749+701 


6, 


.8 


11, 


.8 


8. 


.0 


2, 


.29E+04 


46. 


,79 








42.24 


38 


1758+388 


24, 


.2 


1, 


.0 


1. 


.3 


1, 


.43E+06 


46. 


,73 








42.44 


39 


1800+440 





.9 


28, 


.2 


46. 


.7 


1, 


.91E+03 


45. 


,73 








43.28 


36 


1803+784 


10, 


.6 


1, 


.1 


1. 


.6 


1, 


.04E+05 


46. 


,17 


44, 


.56 


1 


41.49 


25 


1807+698 


9, 


.8 


10. 


.5 


5. 


.0 


6, 


.15E+02 


45. 


,21 


41, 


.40 


12 


40.96 


28 


1823+568 


10 


.4 


4, 


.7 


5. 


.5 


2, 


.95E+04 


46. 


,82 


43, 


.32 


1 


43.28 


25 


1828+487 


3, 


.8 


15, 


.0 


15. 


,2 


1, 


.30E- 


^04 


46. 


,55 


45, 


.26 


1 


44.65 


25 


1830+285 


22, 


.0 


6, 


.9 


1, 


,8 


1, 


.79E- 


f-05 


47. 


,65 


45, 


.41 


1 


42.97 


25 


1845+797 


41, 


.9 


3, 


.8 


0. 


,9 


1, 


.02E- 


f04 


46. 


,33 


42, 


.99 


5 


41.11 


5 


1921-293 


25, 


.5 


9, 


.0 


1, 


,1 


3, 


.90E- 


f05 


48. 


,01 


43, 


.67 


18 


<43.35 


32 


1928+738 


5, 


.3 


11, 


.5 


10, 


,8 


6, 


.88E- 


^03 


46. 


,26 


45, 


.18 


1 


42.02 


25 


2007+776 


5, 


.2 


5, 


.1 


8, 


.4 


2, 


.66E+03 


45. 


,79 








41.28 


28 


2131-021 


9, 


.5 


7, 


.8 


5, 


.8 


7, 


.38E+04 


47. 


,27 


43, 


.66 


19 


42.94 


32 


2134+004 


18, 


.2 


2, 


.1 


2, 


.9 


3, 


.41E+05 


47. 


,64 


46, 


.29 


1 


42.46 


25 


2136+141 


10 


.6 


2, 


.6 


4, 


.1 


2, 


.44E+05 


47. 


,60 


46, 


.02 


1 


41.84 


25 


2144+092 


17, 


.4 


2, 


.1 


3, 


.0 


3, 


.30E- 


f-04 


46. 


,65 








42.97 


30 


2145+067 


3, 


.9 


3, 


.6 


6, 


.6 


8, 


.31E- 


^04 


47. 


,23 


45, 


.79 


1 


42.14 


25 


2200+420 


5 


.8 


7, 


.0 


9, 


.3 


6, 


.27E- 


^02 


45. 


,19 


42, 


.38 


20 


39.84 


28 


2201+315 


5, 


.2 


7, 


.1 


10. 


.0 


8, 


.23E+03 


46. 


,31 


45, 


.46 


1 


43.05 


25 


2223-052 


2, 


.1 


33, 


.1 


26. 


.4 


4, 


.53E+04 


47. 


,11 


45, 


.62 


1 


43.59 


25 
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Table 2. Continued. 



Source 


e 


r 


8 


ni 


log J^kin 


log Lblr 


Refs. 


log Loxt.aGHz 


Refs. 




(degree) 






(cm~^) 


(erg s~^) 


(erg s~i) 




(erg s~^) 




2230+114 


2.7 


27.9 


20.3 


5.51E+04 


47.19 


45.89 


1 


43.33 


25 


2234+282 


20.9 


16.8 


0.9 


2.58E+06 


48.85 


44.52 


1 


42.05 


25 


2243-123 


8.4 


12.6 


5.8 


7.80E+04 


47.32 


45.28 


1 


42.70 


25 


2251+158 


3.8 


32.4 


11.5 


1.67E+05 


47.68 


45.68 


1 


43.84 


25 


2345-167 


20.7 


1.6 


2.3 


3.55E+04 


46.47 


44.38 


5 


42.66 


25 



Column (1): lAU name; Column (2): the viewing angle of jet d; Column (3): the Lorentz factor F; Column 
(4): the Doppler factor 5; Column (5): the normalization factor of electron energy distribution m; Column 
(6): the bulk kinetic power of jet Lkini Column (7): the total luminosity in broad emission lines Lblr, 
also see Chen, Gu & Cao (2009) ; Column (8) : the references for flux of broad emission lines used to estimate 
^blr; Column (9): the radio extended 5 GHz luminosity I/ext,5GHz, from ATCA images; Column (10): 
the references for the radio extended flux. 

References: (1) Cao & Jiang (1999); (2) Scrgccv ct al. (1999); (3) Scarpa & Falomo (1997); (4) Ho ct 
al. (1997); (5) Liu, Jiang & Gu (2006); (6) Wang, Lu & Zhou (1998); (7) Rokaki ct al. (2003); (8) Our 
unpublished measurements of Mg II line; (9) SDSS spectra; (10) Brothcrton (1996); (11) Kuraszkicwicz ct 
al. (2004); (12) C97; (13) Marcha ct al. (1996); (14) Fan, Cao & Gu (2006); (15) Wilkes (1986); (16) Walsh & 
CarswcU (1982); (17) Cao (2000); (18) Jackson & Browne (1991); (19) Rector & Stocke (2001); (20) Corbett 
et al. (1996); (21) Cooper, Lister & Kochanczyk (2007); (22) Wills & Browne (1986); (23) Kharb & Shastri 
(2004); (24) Gelfand et al. (2005); (25) Cao & Jiang (2001); (26) Vermeulen & Cohen (1994); (27) Marshall 
et al. (2005); (28) Perlman & Stocke (1993); (29) Pedlar et al. (1990); (30) Antonucci & Ulvestad (1985); 
(31) Stanghellini et al. (1990); (32) Browne & Murphy (1987); (33) Punsly (1995); (34) Saikia et al. (1990); 
(35) Ulvestad ct al. (1981); (36) Zhang & Fan (2003); (37) Perlman et al. (1994); (38) Wu et al. (2007); (39) 
Tinti et al. (2005). 



